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Abstract 

We present a method for computing the spectrum of black hole radiation of a 
scalar field satisfying a wave equation with high frequency dispersion. The method 
involves a combination of Laplace transform and WKB techniques for finding ap- 
proximate solutions to ordinary differential equations. The modified wave equation 
is obtained by adding a higher order derivative term suppressed by powers of a 
fundamental momentum scale ko to the ordinary wave equation. Depending on the 
sign of this new term, high frequency modes propagate either superluminally or 
subluminally. We show that the resulting spectrum of created particles is thermal 
at the Hawking temperature, and further that the out-state is a thermal state at 
the Hawking temperature, to leading order in ko, for either modification. 
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1 Introduction 



Since Hawking's discovery that black holes radiate a thermal spectrum, |IJ, various other 
derivations of this effect have appeared, || H All seem to depend in some crucial 
way on the very high frequency behavior of the theory. Clearly such a derivation cannot 
be trusted without taking into account backreaction effects from the spacetime metric. 
This has led a number of authors, ||, ||, [7|, ||, to consider the effects of high frequency 
dispersion on the Hawking spectrum. All have found, in the context of specific models, 
that the Hawking radiation remains almost exactly thermal at the Hawking temperature. 

In this paper we show how the leading order contribution to the Hawking flux can 
be obtained by analytical methods for two models containing high frequency dispersion, 
one in which the high frequency modes propagate subluminally and one where the high 
frequency modes propagate superluminally. The method involves a combination of WKB 
and Laplace transform techniques to solve the modified wave equations, similar techniques 
were also used by Brout, Massar, Parentani, and Spindel || for different sub luminal type 
models. We show for both models that the Hawking flux remains exactly thermal at 
the Hawking temperature to leading order in inverse powers of k . We further show, to 
leading order in k , that static observers far outside the black hole see the in- vacuum as 
a thermal state at the Hawking temperature, also in agreement with the ordinary wave 
equation case ||. 

The specific models of high frequency dispersion considered in this paper are obtained 
by adding a higher derivative term, suppressed by a new fundamental momentum cutoff 
ko, to the ordinary wave equation with the appropriate sign to generate either subluminal 
or superluminal propagation of high frequency modes. The low frequency modes however 
behave as in the ordinary wave equation. Various subluminal theories have already been 
considered in ]5|, ||, [7|, ||, in particular the subluminal equation considered in this paper 
is the same as considered in H. Recently Unruh jnj has also considered a superluminal 
modification (although different from the one considered here) to the ordinary wave equa- 
tion. He has shown by numerically solving the modified wave equation that the spectrum 
is very nearly thermal at the Hawking temperature, and that the Hawking particles arise 
from vacuum fluctuations inside the horizon. Since these vacuum fluctuations would in 
principle evolve out of the singularity, a boundary condition at the singularity would be 
required. To avoid this problem, Unruh instead demanded that the vacuum fluctuations 
were in the ground state outside the singularity. In this paper we consider a similar su- 
perluminal modification to the ordinary wave equation, with the same type of boundary 
condition. 

The remainder of this paper is as follows. We begin by introducing the model in 
section ^. In section |^ we discuss the method used to compute the particle creation in 
this model. Sections |] and [5] then describe how the relevant solutions for particle creation 
are obtained for the subluminal and superluminal dispersion relations respectively, and 



I 



in section ^| we end with some conclusions. We use units with c = h = 1. 



2 Model 

We consider a real scalar field propagating in a 2-dimensional black hole spacetime with 
metric 

ds 2 = -dt 2 + (dx - v(x)dt) 2 . (1) 
This is a generalization of the Lemaitre line element of Schwarzschild spacetime where 



v(x) = —y2M/x. We follow the convention that v(x) < and (in units where c = 1) the 
horizon is located at v(xh) = — 1. The action for the field is given by 



S= l -i (Fx 



{{d t +vd x )iP) 2 +iPF(d x )ip]. (2) 



To motivate this action we note that the black hole defines a preferred frame, the frame 
of freely falling observers. In the Lemaitre coordinate system, {d t + v(x)d x ) is the unit 
tangent to free fall observers who start from rest at infinity, and d x is its unit, outward 
pointing normal. Our action comes from modifying the derivative operator only along the 
unit normal d x . 

In the ordinary, minimally coupled action, F(d x ) = d 2 . In this paper we take 



F ± (d x ) = d 2 x ± ±d x . (3) 
k 



The essential difference between these two derivative operators is that F + produces sub- 
luminal propagation of high frequency modes| whereas _F_ produces superluminal prop- 
agation. The simplest way to see this behavior is to look at the respective dispersion 
relations. Varying the action produces the equations of motion 

(d t + d xV )(d t + vd x )<p = d 2 x <\> ± ia^>. (4) 

K o 

Assuming for simplicity that v(x) is constant, we solve ([|) by looking at mode solutions 
of the form 

(pit, x) = exp(i(ut — kx)). (5) 
This produces the dispersion relations 

(w - vkf = k 2 T (k/k )\ (6) 



2 Actually F + produces both subluminal and superluminal propagation, however, the superluminal 
behavior does not play a crucial role in the analysis of this paper. It could in fact be removed by 



considering a model like Unruh's, 0. 
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In figure [I] we plot the square root of (|G[), that is we plot (u — vk) and yk 2 ± k A jk\ as 
functions of k for a fixed lu, along with the square root of the dispersion relation resulting 
from the ordinary wave equation. The intersection points are the allowed wavevector 
roots to (0). The value of the slope of the \jk 2 ± k A jk\ curve evaluated at an intersection 
point is the locally measured (by a freely falling observer) group velocity of a wavepacket 
centered about that wavevector. From figure |l]we see for the \jk 2 + k A jk\ (corresponding 

to F_) curve that this slope is approximately one when k -C k and increases with 
increasing k, therefore the high frequency modes propagate superluminally. Similarly the 
\Jk 2 — k A jk\ (corresponding to F + ) curve has slope approximately one when k <^ ko and 
decreases with increasing k until at some finite k it goes to zero, with a further increase 
in k the magnitude of the slope increases and eventually becomes much larger than one. 
Therefore the "large" wavevector modes propagate subluminally, but the "very large" 
wavevector modes, i.e., with k near ko, propagate superluminally. As already mentioned, 
the latter superluminal behavior is not essential for what we discuss in this paper, it could 
be removed by considering a dispersion relation like Unruh's || in which the slope of the 
dispersion curve asymptotes to zero as k goes to infinity. 




Figure 1: Plot of the dispersion relations for the ordinary wave equation and its subluminal 
and superluminal modifications. The intersection points with the (w — vk) line are the possible 
wavevector roots. For the particular v shown, there are four real wavevector roots for the 
subluminal dispersion relation, but only two real roots for the superluminal case. The other two 
roots in this case are complex. 

An important property of the action (0), when generalized to a complex scalar field, 
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is that it is invariant under constant phase transformations of the field. This leads to a 
conserved current j M . The time component j°, when integrated over a spatial slice, serves 
as a conserved inner product when evaluated on solutions to the equation of motion 
For the metric (|l|), the inner product takes the form 

(F,G)=i Jdx(F*(d t + vd x )G-G(d t + vd x )F*), (7) 

where F(t, x) and G(t, x) are solutions to (f|). Two classes of solutions to the field equation 
(|j) are of interest. The first are the positive free fall frequency wavepackets. They can be 
written as a sum of solutions satisfying 

{d t + vd x )F(t, x) = -iu'F(t, x) (8) 

where u' > 0. The second are the positive Killing frequency wavepackets. These are a sum 
of solutions of the form e~ lu}t f(x) where u > 0. A positive free fall frequency wavepacket 
need not have a positive norm under (0) in general, but does when v(x) = constant. A 
positive Killing frequency wavepacket also need not have positive norm in general, but 
does when v(x) = 0. 

To quantize the field we assume that 4>(x) is a self-adjoint operator solution to the field 
equation which satisfies the canonical commutation relations. We define for a normalized 
positive free fall frequency solution f(t,x) the annihilation operator a(f) by 

a(f) = (fj). (9) 

We make a similar definition for the annihilation and creation operators for a normalized 
positive Killing frequency solution g(t,x) (assuming it is also positive norm). 

3 Computing the Particle Production 

The standard method of computing the amount of particle production in a given wavepacket 
is to propagate this wavepacket back in time to the hypersurface where the quantum state 
boundary condition is defined. In our case we are interested in computing the particle 
production in a late time, outgoing (right-moving), positive Killing frequency wavepacket. 
We assume that the state of the field is the free fall vacuum, which is defined by a(p) |ff) = 
for all positive free fall frequency modes p on the early time hypersurface. Denoting our 
late time wavepacket by ipout, one may show that its number expectation value in the free 
fall vacuum is 

where ip- is the negative free fall frequency part of ip out after being propagated back to 
the early time hypersurface (see || for a detailed derivation). 
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Rather than solving the full equation of motion, (|J), we shall instead restrict ourselves 
to mode solutions of the form 

i/;(t,x) = e~ iwt (j)(x). (11) 

Substituting into (||) (and setting k = 1) produces the ordinary differential equation 
(ODE) 

± (() {iv \x) + (1 - v 2 (x))(p"(x) + 2v(x)(iuj - v'{x))(j)'{x) - iu{iu - v'{x))(j){x) = (12) 

where the ± refers to F± respectively and we have used a prime (/) to denote a derivative 
with respect to x. Restricting ourselves to mode solutions (|TT|) therefore has the advantage 
that we need only solve an ODE. 

To determine the boundary conditions for ([12]) it is necessary to study wavepacket 
propagation in these models. For the subluminal, F + , equation this has been discussed 
in great detail in || ||. The conclusion is that the late time, positive Killing frequency 
packet comes from a pair of ingoing, short wavelength packets, located far outside the 
black hole, and nothingP] else. In particular, no part of the wavepacket piles up against 
the horizon as with the ordinary wave equation, and furthermore nothing comes from 
across the horizon. This led || to conclude that the boundary condition for (|1^) in the 
subluminal case is that the solution decays across (and inside) the horizon^. 

For the superluminal, F_, equation we refer the reader to |K| |TTJ for the details of 
wavepacket propagation. The conclusion is that the late time, positive Killing frequency 
packet comes from a pair of right-moving, short wavelength packets located far inside 
the horizon. This is hardly surprising given that the F_ dispersion relation produces 
superluminal wave propagation. The important point concerning the boundary conditions 
though is that the only packet ever outside the horizon is the late time, positive Killing 
frequency packet^. This leads to the boundary conditions for the superluminal ODE that 
for x ^> the solution reduces to a single mode with wavevector corresponding to the 
outgoing wavepacket. 

Once we have solutions to the ODE's satisfying the boundary conditions just discussed, 
we may easily extract the particle creation. In the subluminal case, the solution at x 3> 

3 There could also be an ingoing, long wavelength packet which arises from scattering as the outgoing, 
positive Killing frequency packet is propagated back in time. For a slowly varying background spacetime 
we expect that the amount of scattering is negligible, which is indeed what we find. 

4 The solution cannot vanish inside the horizon because the coefficient of the highest derivative term 
does not vanish, as it does with the ordinary wave equation. 

5 This is not quite correct. Propagating the outgoing, positive Killing frequency packet back in time, a 
piece of it will scatter into an early time, ingoing packet. We again expect that the amount of scattering is 
very small for a slowly varying background spacetime, and in the leading order approximation discussed 
below we can ignore it. This problem can be dealt with by instead looking for mode solutions which 
correspond to propagating positive free fall frequency packets forward in time, which are necessary for 



computing the full quantum state. This will be discussed in subsection 5.1 
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(where v{x) ~ constant) can be decomposed as 



0(*)=]>>Me^, (13) 
1=1 



where ki(u) are the roots to the subluminal dispersion relation From figure [I] it is easy 
to see that two of these roots are positive and the other two are negative. The late time 
wavepacket corresponds to the small, positive wavevector (k +s ), and the early time, ingo- 
ing wavepackets to the large, positive wavevector (k + ) and the large, negative wavevector 
{kJ) respectively. The small negative wavevector (fc_ s ) corresponds to a long wavelength, 
ingoing wavepacket which will not be important in this leading order calculation. The 
number expectation value for a mode of Killing frequency uj is then 

K(fc-K(Mc 2 _MI (u 

1 ] \u'(k +s )v g (k +s )cUu)\- { ' 

The kinematic factors v g {k) and u'(k) are the group velocity as measured by a static 
observer and the frequency as measured by a freely falling observer of a wavepacket 
narrowly peaked about wavevector k. 

The superluminal equation can be handled in almost exactly the same manner. Far 
outside the horizon, where v(x) is approximately constant and satisfies > v(x) > —1, it 
is again easy to see from figure [I] that there is one positive wavevector root (k +s ) of the 
dispersion relation (|[) (with the plus sign) and one negative wavevector root (k- s ). Our 
boundary conditions dictate that the solution at x ^> is[] 

c+s e ik +^ )x . (15) 

To avoid dealing with the singularity, we shall also assume that v(x) becomes constant 
behind the horizon. When the slope of the straight line in figure [TJ is larger than one, it 
is easy to see that there is one positive wavevector root (k + ) to the dispersion relation 
(§) (with the positive sign) and three negative wavevector roots which we denote as fc_ s , 
fc_ m , and k- in order of increasing magnitude (s denoting small and m denoting middle). 
As we shall see below, only the large positive and large negative wavevector solutions will 
contribute to the solution. It follows that the solution in this region is of the form 

c+e iM")* + C _ e ifc-Hx_ ( 16 ) 



The number expectation value in this case again becomes (14), with the kinematic factors 
appropriate for the superluminal equation, F_. 

6 Note that we use the same notation for denoting the wavevectors for both the subluminal and 
superluminal cases; however, the actual values of the wavevectors for a given u> differ between the two 
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4 Approximate solutions to the subluminal equation 



The methods applied to find approximate solutions to ODE fll2|) are the same as those used 
in solving the Schrodinger equation for a tunneling potential. For the general potential 
one may find approximate solutions to the Schrodinger equation by the WKB method 



14j; however, about a classical turning point (i.e., where the kinetic energy vanishes) the 
WKB approximation breaks down. Approximate solutions can nevertheless be obtained 
by expanding the potential V(x) in the full Schrodinger equation about the classical 
turning point x tp . Solutions to the resulting equation are straightforward to find and are 
valid in the region 

\x - a^r 1 « \n\V\x tp )/V^\x tp )\ (17) 

if V"(x tp ), • ■ ■ , V^ n ~ l \x tp ) all vanish. 

If the WKB solutions are valid in regions on either side of the classical turning point 
which overlap with the region of validity of the solution to the linearized potential equa- 
tion, then we can obtain approximate solutions to the full Schrodinger equation over the 
entire range of x. We now apply this method to ODE (|T2|) for the subluminal equation, 
F + . We begin by finding the approximate solution about the horizon relevant for particle 



creation in subsection |4.1| , and then match this solution to WKB solutions outside the 
horizon in subsection O . From this solution we then compute the amount of particle 
creation in subsection [4.3| . Finally we show in subsection [4.4| that this method can be 
extended to computing the complete out-state to leading order in ko, i.e., compared to 
just computing number expectation values. 

4.1 Approximate solutions about the horizon 

In words the calculations we present in this subsection are as follows. We find approximate 



solutions to ODE (|12|) about the horizon by linearizing the x-dependent coefficients in 
the equation and solving the resulting equation by the method of Laplace transforms. 
These solutions are given by contour integrals in the complex s-plane (where s is the 
Laplace transform variable). The contour Cq (see figure H) corresponding to the x-space 
solution inside the horizon is chosen so that the boundary conditions described in section 
H] are satisfied, i.e., that the solution decays inside the horizon, as these are the boundary 
conditions relevant for particle creation. The x-space solution outside the horizon must 
then arise from a contour that is deformable to Cq. This contour is broken up into three 
separate contours: C\, C 2 , and C 3 (see figure |3|). By comparing the solutions corresponding 
to these contours to the WKB solutions (computed in subsection |4~2|) of ODE (|T2|), we 
show that C3 corresponds to the late time, outgoing Hawking particle and that C\ and 
C*2 correspond to the ingoing, large positive and negative wavevector packets respectively 
from which the outgoing Hawking particle arises. Once this has been done it is a simple 
matter to compute the particle creation as discussed in section 0. 
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We first linearize v(x) and v'(x) about the horizon as 



v(x) ps — 1 + kx (18) 
v'{x) «k + k\x (19) 

where k is the surface gravity of the black hole described by the metric (]!]) and k% is a 
higher order correction to v(x). Substituting into flT2|) and keeping only linear terms in x 
yields 

( ™ ) (x) + 2/tx0 / '(x)+2(-(icj-K) + ^ = 0. (20) 

Validity of this equation requires that \kx\ <C 1 and \k\x/k\ < 1. To leading order we 
may therefore further simplify the equation as 

<p {iv) (x) + 2nx<p"(x) -2(iu- k)4>\x) - iu(iu - k)</>(x) = 0. (21) 

This is the equation we shall use in this paper to find approximate solutions about the 
horizon; however, to compute correction terms to the flux, we must keep at least all linear 
terms in x and possibly even higher order terms in x. We shall discuss this further in the 
conclusions section. 

We use the method of Laplace transforms, [13, I6| , to solve (|2T|). Writing the solution 
clS db Laplace transform, 

<P(x) = [ dse sx 4>(s), (22) 
Jc 

(where C is the contour of integration) and substituting into ( pT|) yields the s-space[] ODE 

d, (M^M)) = s4 " 2( "'^7 M "' K) - (23) 

Equation (p3|) is easily solved as 

*')-^-p(s(t + ^)) <"» 

To obtain the x-space solution, we substitute <j>(s) into (P2] ) and integrate. The choice 
of contour C over which we integrate is dictated by the boundary conditions discussed in 
section i.e., we want 4>(x) to decay inside the horizon. Before finding the appropriate 
contour to produce this behavior let's first understand the generic properties that the con- 
tour must satisfy. Specifically, note that </>(s) is dominated at large \s\ by the exp(s 3 /(6/t)) 
term, and therefore for the integral to converge (assuming the contour C runs to infinity, 

7 We have dropped a boundary term to obtain (^3|). This term will vanish by our choice of contours 
below. 
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which it need not) the contour must asymptote to a region where Re(s 3 ) < since k is 



real and positive. Writing s 
of the three regions 



re 



this implies that the contour must asymptote to any 



7V 7V 

Region 1 <-> — < 9 < — 
6 2 

r> • o 5-7T 7-7T 

Reqionl «-> — < U < — 
y 6 6 

_ 37T 117T 

Reqiona <-> — < < . 

y 2 6 



(25) 



In figure ||] these appear as the unmarked regions. 




Figure 2: Diagram of the steepest descent contour Co- The unmarked regions are directions in 
which the contour must asymptote for the integral to converge. The x 's are singularities of the 
integrand and the wavy line is a branch cut. 



To evaluate the contour integral (|22|) we first consider the x < case. We must choose 
a contour that yields a solution that decays with decreasing x inside the horizon, and 
that asymptotes to any of the three regions (fffi) . To find this solution we approximate 
the contour integral fl2"2"| ) by the method of steepest descents, [ffE] . 

First rewrite the contour integral as 



(f)(x) 



c 



ds g{ 



s e 



(26) 
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where 

g( S ) = S -!-W« (27) 

and 

m-'+US+^Y (28) 

To evaluate (EB) by steepest descents we first locate the saddle points of f(s). These are 



given by the roots of df(s)/ds = 0, | 16| , which in this case are approximated by 



s± « ±V-2kx (29) 

(since we are in a region where cj, k <C |kx|). The contours of steepest descent through 
these saddle points are given by Im(f(s) — f(s±)) = and Re(x(f(s) — f(s±))) < 0. 
Using this one may show that the direction of the steepest descent contours through s + 
and s_ are — tt/2 to n/2 and to 7r respectively It is not hard to show that the steepest 
descent contour, Co, through s + asymptotes to regions 1 and 3 as shown in figure |2|. The 
contour integral in this case is now obtained by the standard formula 



Mz) ~ 9(s + ) l J? e */(* + )+-+ (30) 

V \ x f \ s +)\ 

where a + = 7r/2 if we traverse the contour in the direction indicated in figure ^. To lowest 
order in uo and k this reduces to 

o (x) « -v / 2^(-2/€x)~ 3/4 ^ /(2k) exp (-^V2^|x| 3/2 ) . (31) 

We immediately see that this is exponentially decaying with decreasing x (recall that 
x < 0). The contour through s_ produces an exponentially growing solution, hence our 
desired solution is given by the contour Cq. Finally, note that g(s) is singular at s = 0, 
and that we must choose a branch cut from this point. We choose the branch cut to run 
along the negative real s-axis. 

At this point the reader may be wondering about the validity of the approximations 
made so far. The steepest descents method requires that |x| 3> 1 while validity of the 



approximate ODE (|2l|) requires that \kx\ <C 1. As long as k <C 1, there is always a wide 
range of x values satisfying both conditions, i.e., 1 <C \x\ 1/k. Such k correspond to 
black hole temperatures Tjj 1 (or T# k if we restore ko). For example, if the Planck 
length is also one in these units (i.e., ko = 1 = 1/Zpz), and k is the surface gravity of a 
solar mass black hole, the inequality on x becomes 1 \x\ 10 38 . Clearly there is no 
problem in satisfying this inequality. It is convenient to keep these numbers in mind for 
later approximations. 

Now we turn to evaluating (f>(x) for x > 0. In principle we must evaluate the contour 
integral (|26|) over the same contour as in the x < case, i.e., Cq. However, by Cauchy's 



10 



theorem we may deform the contour (keeping the endpoints fixed) through any region 
in which the integrand is analytic to a new contour, hopefully one where the integral is 
easier to evaluate. In particular we may deform the contour so that it runs through any 
nearby saddle points so that we may again approximate the integral by the method of 
steepest descents. In fact most of the work for these saddle points has already been done. 
They are still given by (|29|) except that now x > and therefore are both imaginary (to 
leading order), i.e., s± ~ ±i\/2K,x. The direction of the steepest descent contours through 
s_ and s + are now given by 7tt/A to 37r/4 and 57r/4 to 7r/4 respectively. From this one 
can easily see that the steepest descent contour, G\ y through s_ asymptotes to regions 2 
and 3 (^) and the steepest descent contour, C 2 , through s + asymptotes to regions 1 and 
2 (|25|), as shown in figure |3|. Evaluating the leading order contributions to these contour 
integrals as before results in 



(j) 2 {x) 



exp 



-i-V2nx 3/2 
3 



exp [ i-y/2~K,x 3 ^ 2 



« e -*T/4 e ™/(2«) (2kx) - 3 / 4 - W( 2k ) 

where we have chosen the directions of the contours as depicted in figure 131. 



(32) 
(33) 




Figure 3: Diagram of the steepest descent contours Ci, C2, and C3. C\ and C2 pass through 
the saddle points s+ and s_ respectively. The unmarked regions are directions in which the 
contour must asymptote for the integral to converge. The x 's are singularities of the integrand 
and the wavy line is a branch cut. 
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The contour C\ + C2 is not by itself deformable to Co, but if we add in the contour 
C3 (see figure ^) which asymptotes to region 2 on either side of the branch cut, then 
C\ + C 2 + C 3 is deformable to C . To evaluate the contour integral over C 3 , first define 
the new integration variable e~ l7T t := sx. This produces 



x 



iio [ K 



C 3 



rft(-t) 



— 1— iuj/n 



exp — t 



— K)X 
2«t 



(34) 



where the new contour C3 runs from infinity just above the postive real axis, counter- 
clockwise about the origin, and back to infinity just below the positive real axis. Ignoring 
the t 3 and t~ x terms in the exponent for the moment, we note that the remainder is just 



a gamma function, i.e., using the integral representation |T7 

1 



i2sin(nv) Jc 3 



dt (-t) 



-l+v. 



we arrive at 



<f) 3 (x) « -2sm\\(TTuu/K)T(-iuu/K)x iu),K . 
To see when this approximation holds, expand 



exp 



J) 



Qkx 3 



+ 



K 2 X 6 



(35) 



(36) 



(37) 



Evaluating the integral 



with the t 3 term produces the correction term to 03 (x) 
smh.(iruj / k) 



3kx 3 



T(-iuo/K)x iullK . 



(3* 



Using the identity T(z + 1) = zT(z) one may show that \5(p^{x) / (p^{x)\ 1 holds if 



1 < \kx 3 



(39) 



If we expand exp(- 
correction term to 



K)x/(2Kt)) in the same way and evaluate the leading order 



-iuo[iuo 

^(x) as before, we find that we need 
\ujx\ <^ 1. 
is well approximated by 



(40) 

(in the asymptotic ex- 



(x) to the mode equation 



When these conditions hold, 
pansion sense). 

To summarize, we have found an approximate solution 
(|T2]) satisfying the boundary conditions that it decay inside the horizon. Just outside the 
horizon the solution is given by <fi(x) = 4>\{x) + 2 (x) + 4>z{x) (see (|32|), (|33|) , and (|36|) 
respectively). We now propagate this solution out to a region where v(x) is essentially 
constant by patching onto WKB solutions which are valid outside the horizon. Knowing 
the solution in the constant v (x) region will then allow us to easily extract the particle 
flux. 



12 



4.2 WKB solutions 



To find approximate solutions to the mode equation (|1"2D by the WKB method, assume a 
solution of the form 

=e ijdxk(x) ( 41 ) 

where the wavevector k(x) is an unknown function of x. Substitution yields 

k 4 - (1 - v 2 )k 2 - 2vuk + u 2 = 

i^-(2k 3 - (1 - v 2 )k - ujv) + Ukk" + 3(k') 2 ) - ik"' (42) 
ax 

where we denote derivatives with respect to x by primes (/). If v(x) is a slowly varying 
function of x, then we expect k(x) also to be slowly varying. We therefore try to solve for 
k perturbatively in derivatives of v(x). More precisely, let x — > ax (we will take a = 1 at 
the end), then (f4*2| ) becomes 

k A - (1 - w 2 )A; 2 - 2vuk + u 2 = 

1 A(2A; 3 - (1 - v 2 )fc - cuw) + \(4kk" + 3(k') 2 ) - —ti". (43) 
a ax a 2 a A 

In this equation we see on the right-hand-side that derivatives of k(x) and v(x) with 
respect to x are suppressed by powers of 1/a. 

Now assume that k(x) may be expanded in inverse powers of a as 

k(x) = k {0) (x) + - k {1) (x) + • • • . (44) 
a 

Substituting into (fE|) and demanding that the coefficients of each power of 1/a separately 
vanish produces an infinite set of equations, the lowest orders being 

(^(o)) 4 _ (l _ v 2 )(k^) 2 - 2vcok + u 2 = (45) 

jfe(i) = 1 A. \n(2(k^) 3 - (1 - v 2 )k^ - vu). (46) 
2 ax 

Although the leading order equation for k^°\ ([45]), can be solved exactly producing a set 
of four wavevectors, the expressions are quite unwieldy. Fortunately, since we are mainly 
interested in Killing frequencies u satisfying w < 1, we only need find approximate 
wavevector roots to this equation. Once these roots are known, the 1/a corrections to 
them can be found by substituting the respective wavevector root k^ into (|46|) and solving 
for k^ 1 '. These computations produce the wavevectors 

k ± = ±VT^+^ + i-^Hl-v 2 )+0(u 2 ) (47) 

Ks = ^- + 0(u 3 ) (48) 
1 + v 

k- s = --^ + 0(cu 3 ), (49) 
1 — v 
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where we have set a = 1. The corresponding WKB solutions are 



(j) ± (x) ps (1 — v(x) 2 )~ 3 ^e ±i I dx V 1 '^ 2 e iuJ I dxv ( x )/( l - v2 i x )) 

<P +S {X) ~ e iufdx/(l+v(x)) 

~ e -iufdx/(l-v(x)) m 



(50) 
(51) 
(52) 



k { ±\x) 


3 


w(x)f '(x) 


k ( ±\x) 


~ 2 


(1 - V 2(a;))3/2 



The condition of validity for these approximate solutions is that \k^\x) / k^ {x)\ <C 1. 
For the fc± wavevectors this ratio is 



(53) 



We are interested in v (x)'s containing black holes. The horizon of a black hole in these 
units is located at v(xh) = — 1, therefore the right-hand-side of fl53D clearly becomes 
arbitrarily large as we approach the horizon (assuming that v'(x) ^ which are the 
only cases we consider here). It follows that the WKB approximation will break down 
around the horizon. Far from the horizon (and outside the black hole) v(x) asymptotes 
to a constant — 1 < vq < and v'(x) goes to zero, therefore the WKB approximation 
will be valid. For the k +s mode a ratio similar to ( |53"|) holds, and therefore the WKB 
approximation again fails around the horizon, but is valid far outside of it. To compute 
this ratio, we must compute k +s to order 0(u 3 ), however since we will not need this later 
we do not give the explicit expressions here. 



4.3 The spectrum 

We now have all the ingredients necessary to compute the leading order spectrum of black 
hole radiation. First note that validity of the WKB approximation (|53|) requires that 
1 V kx 3 . Furthermore the approximate solution from the Laplace transform method is 
valid when 1 <C \x\ <C 1/k. Since we are considering cases where k<1, there is always 
a region where both the WKB and Laplace transform solutions are valid. 

If we evaluate the integrals appearing in the WKB solutions of (|50"D, fl5"T|) and, (|52| ) 
respectively in a region just outside the horizon where v(x) is given by the linearized 
expression (|T8D, we find that the solution outside the horizon obtained from the Laplace 
transform method can be expressed as 

(f){ x ) = e - i7r/4 v / 2^(-ie-^ /(2K) 0_(x) + e™ /(2K) 0+(x)) 

- e -™ /K (e 2 ™ /K -l)Y{-iu/n)(j) +s {x). (54) 

With <j)(x) decomposed in terms of the WKB solutions, we are allowed to evaluate it at 
large x, i.e., where v(x) is essentially constant. In this region the WKB solutions reduce 



14 



to simple modes (up to multiplicative constants) and we need only extract the coefficients 
of these modes in order to compute the particle creation rate as given in (|14]) . A simple 
computation yields 

N ( u ) = e w*-i ' w 

exactly a thermal spectrum at the Hawking temperature T# = k/(2tt). 

As a check on our results, recall that invariance of the action (|2[) under constant phase 
transformations of the field leads to the existence of a conserved current j M . When this 
current is evaluated on fixed Killing frequency mode solutions, the time component j l 
is manifestly time independent and the conservation law reduces to d x j x = 0, i.e., j x is 
constant. The exact form of j x is complicated, but when evaluated on a mode of the form 
c(u) exp(— zcut + ik{u)x) in a region where v(x) is constant, it reduces to 

f{k{u)) = uj'(k(u]))Vg(k(uj))\c(k(u)))\ 2 . (56) 

For the solution of interest to us, i.e., the one that decays inside the horizon, the spatial 
part of the current must vanish everywhere. Using the solution given by ([54]), it is easy 
to show that this is indeed the case. 

We have made a number of approximations in order to compute the leading order flux 
given by fl55|), we would now like to collect them to find out what restrictions they place 
on the allowed parameter ranges. First recall that we have restricted k to values k C 1 
(in units of ko — 1). Physically this says that we only expect to get a thermal spectrum 
of radiation when the black hole is large, and therefore has a small temperature. 

The range of validity of the Laplace transform solutions in the spatial variable x is 
also restricted. We have already seen that we need \x\ ^> 1 for the steepest descents 
approximation to hold for the various contour integrals and \kx\ <C 1 for the approximate 
ODE fl2Ip to be valid. Closer investigation of the correction terms to the WKB solutions 
and the Laplace transform solutions shows that we need 

l//t 1/3 < \x\ < 1/k 3/5 . (57) 

This inequality means that the matching of the WKB and Laplace transform solutions 
can be done anywhere in this range. To derive these inequalities, note that the WKB 
and Laplace transform solutions cannot be matched to arbitrary order because in fact 
they solve different equations. The Laplace transform solutions were obtained by finding 
approximate solutions to the linearized v(x) ODE (pT|), while the WKB solutions were 
obtained by finding approximate solutions to the full ODE. By computing the first set of 
correction terms to both the WKB and Laplace transform solutions that differ, i.e., do 
not match, and demanding that they are small (compared to the leading order terms), 
we derive the above inequality. 

These correction terms also restrict the range of uo. Since the corrections are x de- 
pendent, by matching the WKB and Laplace transform solutions about an appropriate 
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x satisfying (|5T ) we can minimize the difference between the solutions. To carry this 
out we compute the leading order correction terms for the WKB and Laplace transform 
solutions that differ and sum their absolute values. As a specific case, the leading order 
relative difference between the WKB solution corresponding to the wavevector k + and 
the Laplace transform solution corresponding to the contour C2 is 

ff ~U 2 J K (2/tx) 3 / 4 + 40 k { } 

in the limit of u 3> k. Minimizing this with respect to x we find that K,x m i n ~ (9klu 2 / (47r)) 2 / 13 

Demanding that a Cl we arrive at the constraint 

I nm 13/10 

A similar computation can also be carried out for the WKB and Laplace transform solu- 
tions corresponding to the and k +s wavevector roots. In the limit of u ^> k a weaker 
constraint than above is obtained from the k +s root, and the same constraint as above 
is obtained for the wavevector. In the same manner as above a lower bound on the 
allowed range of ui can be obtained. From the k + wavevector we find that we need 

It follows that the range of validity of the analytical results presented above are bounded 
both above and below in the parameter to. 




4.4 Computing the quantum state 

We have so far computed the outgoing flux of particles from a black hole for the sublumi- 
nal equation of motion (|4]) (with the plus sign). There is, however, much more information 
contained in the quantum state than just number expectation values, for instance, cor- 
relations. It is therefore of interest to compute the full quantum state in this modified 
theory and compare it to the state arising with the ordinary wave equation, as already 
computed by Wald 0. 

To be a bit more precise, we are not actually "computing" the quantum state because 
we already know what it is, i.e., we have assumed that it is the free fall vacuum. What 
we are going to do is re-express this state in terms of a vacuum state defined by late 
time observers. We define the out-Hilbert space as the tensor product of Hilbert spaces 
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on either side of the horizon. Outside the horizon we use Killing frequency to define the 
Hilbert space, as we have done so far. Inside the horizon we do the following. We take 
a v(x) that asymptotes to a constant (less than —1), then we define our Hilbert space in 
this region using free fall frequency. If we are only interested in the observations made 
by the outside observer then we would trace over the inside degrees of freedom, in which 
case the Hilbert space we use inside the horizon is irrelevant. For the computations here 
though, the choice we have made is the simplest. 

Our method of computing the state is very similar to the techniques employed by Wald 
0. In the time dependent picture we would take an ingoing, positive free fall frequency 
wavepacket and evolve it from the hypersurface where the free fall vacuum is defined to 
the hypersurface where the out vacuum is defined. If this packet is sufficiently peaked in 
its wavevector, we may follow it on the dispersion relation as discussed in detail in H. It 
is not hard to see that this packet will propagate toward the horizon and scatter (mode 
convert), the reflected piece propagates (forward in time) away from the horizon to the 
region where v(x) is essentially constant, and the transmitted packet propagates deeper 
inside the horizon to where v(x) is essentially constant. The late time packet outside the 
horizon corresponds to the Hawking particle, and the late time packet inside the horizon 
corresponds to the partner. If the initially positive free fall frequency packet ^+fr contains 
only positive Killing frequencies, the final packet ip+out outside the horizon will contain 
only positive Killing frequencies and the final packet ip-m inside the horizon will contain 
only negative free fall frequencies. The annihilation operator associated with i/j+s, i.e., 
a(ip + ff) := (ip + f[,(f)), annihilates the free fall vacuum. Using the time independence of the 
inner product, we therefore derive the equation 

(a(^ +ou t)-a t (^ in ))|ff) = 0. (62) 

Similarly, if the initially positive free fall frequency packet contains only negative 
Killing frequencies, the final packet outside the horizon will contain only negative Killing 
frequencies and the final packet inside the horizon will contain only positive free fall 
frequencies. Using this, a relation similar to (|62| ) can be derived. As shown by Wald ||, 
given a complete set of relations like fl62j) (constructed by taking a complete set of ingoing, 
positive free fall frequency packets), we can re-express the free fall vacuum in terms of 
the out vacuum. 

As we have done thus far, we shall actually use mode solutions instead of wavepackets. 
We derive the mode solutions which, when appropriately summed, produce the time de- 
pendent wavepacket solutions just discussed. We have already derived one mode solution 
in (|54]), although it is not of the form that we want. Rather it decays inside the horizon 
and is a superposition of plane waves with wavevectors k +s , k + , and k_ far outside the 
horizon (where v(x) is essentially constant). In the time dependent picture it corresponds 
to propagating a pair of ingoing, positive and negative free fall frequency wavepackets 
forward in time, with just the right relative weights so that the entire packet completely 
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mode converts around the horizon, turns around, and propagates out to the constant 
v(x) region. To obtain the mode solutions that we want, we shall construct another 
mode solution below which is a superposition of plane waves with wavevectors k + and /c_ 
far outside the horizon (where v(x) is essentially constant), whereas inside the horizon 
(where v(x) is again essentially constant) it reduces to a plane wave with wavevector 
k-. In the time dependent picture this mode solution corresponds to propagating a pair 
of ingoing, positive and negative free fall frequency wavepackets forward in time, with 
just the right relative weights (although different than above) so that the entire packet 
propagates across the horizon and converts into just a negative free fall frequency packet. 
By adding these two mode solutions, call them the n-modes since they are the relevant 
ones for computing number expectation values, with the correct relative coefficient, we 
can eliminate the fc_ (k+) mode outside the horizon. These are the mode solutions we 
want, call them the s-modes since they are the relevant ones for computing the state, 
because they correspond to propagating an ingoing, positive (negative) free fall frequency 
packet forward in time which splits around the horizon into a pair of wavepackets, one 
propagates back away from the horizon, and the other falls inside the black hole. 

Because the details of computing the other n-mode solution discussed above are es- 
sentially the same as discussed in subsections [O] and |4^, we shall only sketch the com- 
putation. We first solve the mode equation (|TJ) (with the plus sign) about the horizon 
by the method of Laplace transforms exactly as before. The only difference is that the 
contour of integration must be changed so as to satisfy the boundary conditions that 
the solution reduce to a plane wave with wavevector fc_ inside the horizon where v(x) is 
essentially constant. A straightforward computation shows that the contour C4 shown in 
figure f| does the job, i.e., evaluating the contour integral (|26| ) over C4 and propagating 
the solution deeper inside the black hole by the WKB approximation to the constant v(x) 
region shows that the solution is 

4 (x) « 2e™ /K sinh(vrc; / '«)r(-iw / V)0_(x), (63) 

which is the boundary condition we want. 

To evaluate the solution outside the horizon, we deform the contour C4 into C 2 and 
C5 as shown in figure (||). These we can evaluate by the method of steepest descents. In 
fact C 2 is exactly the same as the contour C 2 before, see figure C 5 is the same as C\ 
before (see figure 0), except that it lies on a different Riemann sheet, so only the overall 
scale changes. The complete solution outside the horizon after being propagated out to 
the constant v(x) region by the WKB approximation is 

02 (x) + 05 (x) « e- i7r/4 e™ /{2K) (0+ (x) - *e™ /K 0_ (x) ) . (64) 

Combining this with 04 (x) produces the connection formula 

2e™ /K sinh(7rcj / 'k)T{-iuj / 're)0_(:r) <-> e^" /4 e nul/i2ft) (0+(x) - ie nuj/K <f)^{x)) (65) 
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Figure 4: Diagram of the steepest descent contours C2, C4, and C5. C2 and C5 pass through the 
saddle points s+ and s_ respectively, the solutions for these contours are valid for x > 0. The 
solution corresponding to the contour C4 is valid for x < 0. The unmarked regions are directions 
in which the contour must asymptote for the integral to converge. The x's are singularities of 
the integrand and the wavy line is a branch cut. 

where the left- hand- side refers to x < and the right-hand-side to 1 > 0. This is the 
second of the n-mode solutions, the first is given in . 

To obtain the s-mode solutions, we add the n-mode solutions given in fl54|) and (|65|) 
with the correct relative coefficient so that either the 4>-(x) mode or the <p+(x) mode 
cancels at large positive x, this produces the connection formulae 

c *W^_ (a .) «-> +g (x) + e-^ 4 <Mx) (66) 

1 [—lUJ/K) 

e-^y^x) «-> +5 (x) -e— /4 e -W(2i») v^F« 

1 {—tU/K) 

where again the left-hand-side refers to negative x and the right-hand-side to positive x. 
Noting that the annihilation operator associated with the modes 4>+(x) and <p*_(x) at large 
positive x (in the time dependent picture these would be the early time ingoing, positive 
free fall frequency wavepackets) both annihilate the free fall vacuum |ff), we derive the 
following relations as in (^3), 

(at(0 +s )-e^ K a(0*_))|ff) = O (68) 
( a t(0^)_ e -A a(0+s)) | ff) = o . (69) 
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These two relations completely determine the +s and (jf_ content of the free fall vacuum 
for the Killing frequency uj. Using these relations it is simple to show that the free fall 
vacuum is a thermal state at the Hawking temperature, exactly as with the ordinary wave 
equation, 0. 



5 Approximate solutions to the superluminal equa- 
tion 



Since the calculations involved in solving the mode equation ( |12|) with the superluminal 
operator F_ are virtually identical to those given above for the sub luminal operator F + , 
we shall be brief. Computing the approximate solutions outside the horizon by the WKB 
approximation proceeds exactly as before. The main difference compared to the sublu- 
minal case is that we are now interested in WKB solutions both inside and outside the 
horizon. Outside the horizon, the relevant^ WKB mode is 

<j) +s (x)^e iuJ I dx/il+v ^\ (70) 

and inside the relevant WKB modes are 

<f>±(x) « (-1 + v 2 (x) ) ' 3/4 e ±{ f dx y/- 1+v2 W e <*> I ** ^/(i-* 2 (*)) _ ( 71 ) 

As before, it is straightforward to show that these approximate solutions break down 
around the horizon, but far enough outside the horizon they are valid. To find the 
appropriate connection formula for these solutions, we now find an approximate solution 
across the horizon. 

Solutions about the horizon can be obtained by again linearizing v(x) as in flT8|) and 
solving the resulting approximate mode equation by the method of Laplace transforms. 



The linearized equation is just (21) with a minus sign inserted before the fourth derivative 



term. Writing the solution as a Laplace transform as in (|22| ) and substituting into the 
equation produces the s— space ODE (p3| ) with s 4 — > —s 4 . This equation is again trivial 
to solve. One finds upon writing the solution in the form 



(f)(x) = I ds g(s)e xf{s) (72) 
J c 

that 



ZKX \ 6 s 



8 There are other linearly independent solutions both inside and outside the horizon, but they will not 
be needed in this calculation. 
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and 

g (s) = s" 1 -*"/*. (74) 

Before evaluating ([72]), first note that at large |s|, the integral is dominated by 
exp(— s 3 /(6k)), and therefore for the integral to converge, the contour must asymptote to 
a region where Re(s 3 ) > 0. Writing s = re l9 , these regions are given by 

— 7T 7T 

Region 1 <-> — < < - (75) 
ite#ion 2 - < < — (76) 

TlT 37T 

Region?, «-> — < < — , (77) 
6 2 

and are the unmarked regions in figure [|. 

To evaluate the contour integral ([72]), recall from section |3| that our boundary condi- 
tions are specified outside the horizon and state that the solution must reduce to a plane 
wave with wavevector k +s in the constant v(x) region. This solution does not correspond 
to a saddle point because those solutions are either exponentially growing or decaying. 
It is not hard to guess what contour we need though, given our past experience with the 
subluminal dispersion relation. If we take a contour C$ that encircles the branch cut and 
asymptotes to regions 2 and 3 (|7"7|) , see figure §, we get a contour very similar to the 
contour C3 in figure |3|. The approximations that went into evaluating that contour also 
work here. The result is exactly — 3 (x) of (0). This solution is just 0+ s (x) (|70D up to 
a multiplicative constant, and therefore the contour C§ produces the correct boundary 
condition at x 3> 0. 

To evaluate ( [72] ) for x <C we use the steepest descents approximation. The saddle 
points in the integrand of 072]) are given by 



s± « ±is]2K\x\, (78) 

and the steepest descent contours must pass through these points in the directions — 7r/4 
to 37r/4 for s + and 7r/4 to 57r/4 for s_. The contours, C7 and C%, therefore must asymptote 
to regions 1 and 2 and regions 2 and 3 respectively, see figure || Furthermore, Cj + C% 
is deformable to C 6 , and therefore is the contour we want. Evaluating ([72]) over C 7 + C 8 
by the steepest descents approximation and expressing the result in terms of the WKB 
solutions ([Fl]) results in 

7 + <f) 8 {x) w e"^ /4 v / 2^(e-™ /2K 0_(x) - ie™ /2K 0+(x)). (79) 

We now have the complete solution for all x, which can be displayed as the connection 
formula 

e-^ 4 v / 2^(e—/ 2 >_(x) - ie^cj> + {x)) <- ~^ +5 (x) (80) 

1 (1 + zcu/fi;) 
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Figure 5: Diagram of the steepest descent contour C%. The unmarked regions are directions in 
which the contour must asymptote for the integral to converge. The x 's are singularities of the 
integrand and the wavy line is a branch cut. 

where the left-hand-side refers to x < and the right- hand- side to x > 0. Evaluating (|80|) 
at x 3> and i<0 allows us to pull off the coefficients c +s (u) and c_(a>). Substituting 
into fll4]) again produces a thermal spectrum at the Hawking temperature. 

As before we have made a number of approximations to arrive at this result. Collecting 
these approximations together, we can compute the range of validity of these results. 
However, because the difference between the solutions for the subluminal equation and 
the superluminal equation is only the change of a few signs in the end, then the difference 
between the WKB and Laplace transform solutions for the superluminal case is essentially 
the same as the difference between the WKB and Laplace transform solutions for the 
subluminal case. Therefore the constraints on the range of validity of the solutions in the 
parameters k, x, and u are the same as in the subluminal case. 

5.1 Computing the quantum state 

For the superluminal equation we can as well compute the decomposition of the free fall 
vacuum in terms of particle states as seen by late time observers. We define the out 
Hilbert space as before, i.e., we take it to be a tensor product of Hilbert spaces inside and 
outside the horizon respectively. Outside the horizon we define the Hilbert space using 
Killing frequency, and inside we use free fall frequency (as before we take a v(x) that 
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Figure 6: Diagram of the steepest descent contours, C-j and Cg, through the saddle points s + 
and S- respectively. The unmarked regions are directions in which the contour must asymptote 
for the integral to converge. The x's are singularities of the integrand and the wavy line is a 
branch cut. 

asymptotes to a constant smaller than —1 behind the horizon). 

To compute the decomposition of the free fall vacuum, we again look for mode so- 
lutions which when summed together produce an early time positive free fall frequency 
packet propagating toward the horizon (but now located behind the horizon). Around 
the horizon this packet scatters (mode converts) into a pair of packets, a reflected packet 
which propagates deeper inside the black hole to the constant v(x) region and a trans- 
mitted packet which propagates across the horizon out to the constant v(x) region. In 
this picture the transmitted packet corresponds to the Hawking particle and the reflected 
packet to the partner. From such a solution we could obtain an equation for the free fall 
vacuum analogous to (p2]). 

To compute the s-modes (those needed to compute the state) we again first compute 
the n-modes (those needed to compute number expectation values). We have already 
computed one n-mode given by (^). Recall that in the time dependent picture this 
corresponds to a pair of positive and negative free fall frequency packets propagating 
toward to the horizon (and located inside the black hole) with just the right relative 
coefficient that the entire packet propagates across the horizon out to the constant v (x) 
region. The other n-mode therefore correponds in the time dependent picture to a pair 
of positive and negative free fall frequency packets propagating toward the horizon with 
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the right relative coefficient such that the entire packet is reflected and propagates deep 
inside the black hole to where v(x) is constant. The appropriate mode solution therefore 
must decay outside the horizon. 

Computing this mode solution involves the same techniques used already many times, 
so we shall simply quote the result. The connection formula expressed in terms of the 
WKB solutions is 

(e- m/4 v / 2^(-ie m " /{2K) 0+(^) + e 3 ™ /(2K) (M^)) 
+% e ™ /K Wi T ■ / ^ -m(*)) ~ e-^/ 4 e™/ K + (x) (81) 

1 (1 + IUJ I K) 

where the left-hand-side refers to x < and the right-hand-side to x > 0. The solution 
<f)+(x) for positive x decays exponentially with increasing x, and therefore satisfies the 
boundary conditions. 

To compute the s-modes we take linear combinations of the n-modes ( |30D and (|8*I|) 
to eliminate either <f>+(x) or behind the horizon. This results in the connection 

formulae 

+ iNe™ /K 4>_ m {x) «-> iN(j) +s (x) (82) 
0+(x) + iV0_ m (x) <-> iVe™ /K +s (x) (83) 

where 

Ar = e ^ /4 e-^ /(2K) ^ . (84) 

v27tk sinh(7ro;//«) T(l + iuj/k) 

These relations are enough to carry out the decomposition of the free fall vacuum as 
discussed in subsection fO[ In particular the method of obtaining the equations ( j69| ) on 



the free fall vacuum for the subluminal dispersion relation follows exactly in this case as 
well, with the replacement of 0_ in the subluminal case by 0_ m in the superluminal case. 



6 Conclusions 

We have considered two different modifications of the wave equation in a black hole 
spacetime, one producing subluminal propagation of high frequency modes and the other 
superluminal propagation of high frequency modes. We have shown that both equations 
give rise to exactly a thermal spectrum of radiation from a black hole to leading order in 
an expansion in powers of 1/fco- It is natural to try to push the analysis further to obtain 
a correction term to the outgoing flux. We immediately run into the following difficulty 
though. In obtaining an approximate solution to flT2"|) about the horizon, we actually 
solved instead just the linearized equations (|2T|). To obtain a better approximate solution 
about the horizon, we need a better approximation to (0). We could, for instance, keep 
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higher order terms in x when expanding v(x) and v'(x). If we try to solve the resulting 
equation by the Laplace transform method, we find that we get a higher order ordinary 
differential equation in the Laplace transform variable s. In other words, the s-space 
equation is really not much better than the original x-space equation. 

An important assumption made in deriving the thermal radiation for the superluminal 
equation of motion was that positive free fall frequency modes, located behind the horizon, 
were in their ground state. Clearly we don't know a priori whether this is the physically 
correct quantum state condition. In principle we would have to begin with a quantum 
state which evolves into a black hole, and then ask if these modes actually are in their 
ground state. This requires quantum gravity. A more realistic problem to tackle at this 
time is simply to ask where these modes came from in a semiclassical approximation. 
One would guess naively from the singularity. Recent investigations, JF|, |T^] , have shown; 
however, that for certain charged black holes it is possible that these modes simply reflect 
outside the singularity and become ingoing modes, backward in time. This would have 
important implications because it would mean that the Hawking radiation, even for an 
eternal black hole, would originate from ingoing modes, and therefore we would not have 
the infinite degrees of freedom problem, [[?]]. 

We end by noting that the subluminal model (and possibly the superluminal model 
as well) considered in this paper suffers from the "stationarity puzzle". If we try to 
propagate the outgoing modes backward in time all the way out to infinity, where v(x) 
goes to zero, then there can be no particle creation by conservation of Killing frequency. 
One way out of this problem is to introduce time dependence into the equation of motion 
(perhaps via backreaction) to destroy the Killing symmetry. A step in this direction is to 
put the ordinary wave equation on a spatial lattice, this has the advantage of introducing 
naturally a short distance cutoff and at the same time destroying the Killing symmetry (for 
discretizations of most spatial coordinates). Such a model is currently being investigated 
T8fl by techniques similar to those described in this paper. 
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